Phase diagram and magnons in quasi-one-dimensional dipolar antiferromagnets 
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We investigate antiferromagnetic spin chains, which are 
coupled by a weak antiferromagnetic exchange interaction. 
The spins are located on a hexagonal lattice, i.e. frustration 
is present when three-dimensional order sets in. Typical real- 
izations of such systems are the halides ABX3. In this work 
we particularly study the role of the long-range dipolar in- 
teraction within the framework of a Heisenberg model with 
nearest-neighbor exchange and additional dipolar interaction. 
We perform a classical ground-state analysis and show that 
the spin configuration is sensitively dependent on k' , the ratio 
of the dipolar interaction to the interchain interaction, as a 
consequence of their competing character. Several commen- 
surate and incommensurate phases arise in the different re- 
gions of the parameter space. The ground-state investigations 
are supplemented by a stability analysis by means of a linear 
spin-wave calculation. From the magnon spectra we can show 
that all commensurate phases are stable against fluctuations. 
In comparison with experiments (CsMnBr3, RbMnBra) we 
obtain good agreement for the energy gaps. From this we 
conclude that the dipolar interaction is the most important 
source of anisotropy in these Mn-compounds. 



I. INTRODUCTION 

Unconventional magnetic systems have attracted the in- 
terest of experimental and theoretical physicists in the 
last few years [0. In these systems competing inter- 
actions and/or geometric frustration due to the under- 
lying lattice can lead to unconventional ground states, 
magnon spectra and magnetic phase diagrams Fur- 
thermore, fluctuations are enhanced in systems with frus- 
trated ground states as well as in low dimensions. 

An interaction, which is often competing with respect 
to the exchange interaction is the dipole-dipole interac- 
tion In real systems the dipole-dipole interaction 
(DDI) is always present in addition to the short-range 
exchange interaction. Although its energy is lower than 
the exchange energy it plays an important role in low- 
dimensional systems due to its anisotropic and long-range 
character. 

The most famous quasi-one-dimensional systems are 
the ternary compounds ABX3 (A alkaline, B transition 
metal and X halogen), which have been studied inten- 
sively theoretically and experimentally in the context of 
Haldane's phase and solitonic excitations ||. In these 
systems the carrier of the magnetic moment, the B-ions, 



are located on a hexagonal lattice ||. In this work, the 
Mn compounds arc of particular interest. Because the 
angular momentum L is zero, no crystal-field splitting 
occurs in these systems and the dipole-dipole interaction 
should be the most important anisotropy. 

The influence of the DDI in quasi-one-dimensional, an- 
tiferromagnetic spin chain systems has not yet been stud- 
ied very thoroughly. Instead, the DDI is often replaced 
by a single-ion anisotropy JtJ or the coupling between 
spin chains is neglected M. 



II. MODEL 



The Hamiltonian of the dipolar antiferromagnet reads 



(1) 



with spins Si at hexagonal lattice sites x;. The first term 
describes the exchange interaction Jiy which includes the 
intrachain as well as the interchain interaction. In the 
following we consider only nearest-neighbor exchange, i.e. 



Jw 



— J 1,1' along the chains 

— J' 1,1' within the basal plane 



For a hexagonal lattice the Fourier transform of the ex- 
change energy is given by 

J q = —2Jcosq z — 2J' (cosq x + 2cos(-|-) cos(^^q a )^ . 

(2) 

Here and in the following we measure wave vectors in 
chain direction in units of 1/c, and wave vectors within 
the basal planes in units of I /a. The second term in 
Eq. (nl) is the classical dipole-dipole interaction 



ia/3 1/ n2 



5 a 







3(x z - xp 



t (x/ -X// 



X; - X;' 



(3) 



This term is evaluated by means of the Ewald summation 
technique P,M, which allows the consideration of the 
long-range nature of the three-dimensional DDI in terms 
of fast convergent sums. 
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III. GROUND STATES 

In this section we calculate the classical ground states of 
the system as a function of the ratio of the dipolar energy 
to the interchain exchange interaction 



V Z J' 



(4) 



where V z = ^a 2 c is the volume of the primitive cell. 
Due to the quasi-one-dimensionality of the systems (the 
ratio J' /J is 10~ 2 . . . 10~ 3 ) we consider only antiferro- 
magnetic spin configurations along the chain axis, i.e., 
we restrict q z to tt. The ground-state energy reads 



with 



E, 



= - £ SlMtt'Sv SqMjS-q 



l^l< 




(5) 



(0) 



where A q 3 



A q 3 = for q z — tt. The ground states 
are specified by those wave vectors which belong to the 
largest eigenvalues of the matrix M q in Eq. (||). We ob- 
tain the following phases for decreasing n' (for the lattice 
constants we took the values of CsMnBr3 jnj): 

1. Ferromagnetic phase: k' > k{ — 200.50 

In the region where the dipolar energy is large com- 
pared to the interchain exchange, the minimum of 
the ground-state energy is reached at qi = (0, 0, tt). 
This means that spins within basal planes are or- 
dered ferromagnetically, but still antiferromagneti- 
cally along the chains. 

2. Incommensurate phase I: 200.12 = k' 2 < k' < k[ 

The ground state is an incommensurate phase in 
this parameter region. The wave vector moves 
continuously from qi to the wave vector q2 = 
(0, -^i 7 !") (thick solid line in Fig. |]) or any other 
path rotated by 60° (dashed lines in Fig. |]). 

3. Collinear phase: 17.25 — k' 3 < n' < n' 2 

Spins within basal planes are oriented ferromagnet- 
ically in chains, that are aligned antiferromagnet- 
ically to one another (see Fig. |2|). Because of the 
six-fold symmetry, there are six such ground states 
resulting from rotation of the ferromagnetic chains. 
Note that the continuous degeneracy is lifted. 

4. Incommensurate phase II: < n' < k' 3 

In this incommensurate phase the wave vector mo- 
ves from q2 to q3 = (^p,-^=,7r). The incommen- 
surability appears because of the finite slope of the 
dipolar tensor at q3 jl2| and the parabolic behavior 
of the exchange energy. 



5. 120°-structure: «/ = 

This ground state is characterized by a three- 
sublattice spin configuration in each basal plane. 

The coupling of the spin space to the real space induced 
by the DDI forces the spins to align within the lattice 
basal planes for all four phases, in which the DDI is non- 
zero. Thus, the DDI leads to an in-plane anisotropy. 




Q4 = (f,0) 



FIG. 1. Paths of the wave vector that minimize the 
ground-state energy in the Brillouin zone of the hexagonal 
lattice (q z = tt). 

The phase diagram for the whole parameter region of k' 
is shown in Fig. |^. 
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FIG. 2. Spin configurations within the basal plane for dif- 
ferent ratios of dipolar to interchain interaction Only for 
«' = the 120°-structure is established. For infinitesimal «' 
the phase IC II is favored. The spin configurations of incom- 
mensurate phases are not sketched. 



In summary, we find three commensurate and two incom- 
mensurate phases for arbitrary value of k 1 . 



IV. MAGNON SPECTRA 

The spin-wave calculation is of interest in its own right 
and also serves to scrutinize the stability of the phases 
found against fluctuations. To that end we write the 
Hamiltonian (Q) in terms of creation and annihilation 
operators employing the Holstein-Primakoff transforma- 
tion [ p^[ . For low temperatures an expansion up to bilin- 
ear terms can be used, leading to linear spin-wave theory. 
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The Hamiltonian then is diagonalized with a Bogoliubov 
transformation from which we obtain the spin-wave fre- 
quencies. This investigation is restricted to the commen- 
surate phases, because the infinite primitive cell of the 
incommensurate phases resists such an analysis. 
Ferromagnetic phase 

The magnon spectrum of the ferromagnetic phase 
is stable for the parameter region given in Sec. [I]. 



The rotational invariance of the spins around the chain 
axis leads to a Goldstone mode in the spectrum p[ . 
Collinear phase 



The collinear phase of Sec. Ill is also stable against 
fluctuations. There is no Goldstone mode due to the 
discrete degeneracy of the ground state. 
120°-structure 



We argued in Sec. [II that the 120°-structure is unsta- 
ble for infinitesimal dipolar energy due to a linear slope of 
the dipolar tensor at the ordering wave vector. However, 
we performed a spin-wave calculation based on a com- 
mensurate 120°-structure, where the spins are located 
within the basal planes of the lattice |[5) . It turns out 
that the spin-wave spectrum is stable for weak dipolar 
energies, from which we conclude that this commensu- 
rate ground state is a good approximation. 



Thus, we do not need any single- ion anisotropy to ex- 
plain these results. 



V. SUMMARY 

We found three commensurate and two incommensurate 
phases for different values of the ratio of dipolar to inter- 
chain interaction due to the competing character of those 
two interactions. We showed via linear spin-wave theory 
that all commensurate phases are stable against fluctu- 
ations and that the incommensurate phase IC II can be 
approximated by a 120°-structure for weak dipolar ener- 
gies. 

The spin-wave gaps of CsMnBr3 and RbMnBr3 are 
in good agreement with the experiment, which shows 
that the dipolar energy is the most important source of 
anisotropy in these Mn-compounds. 

This work has been supported by the BMBF under 
contract number 03-SC5-TUM and the DFG under con- 
tract number PI 337/1-2. 
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FIG. 3. Magnon spectrum of the 120° -structure perpendic- 
ular to the chain direction. 

The spin-wave frequencies resulting for the dipolar 120°- 
structure for CsMnBr 3 are shown in Fig. [?], where we 
used J = 215 GHz and J' = 0.41 GHz respectively [||. 
This leads to k' = 0.774, i.e., this substance is in the IC 
II region of the phase diagram in Fig. |[ Including the 
DDI, from the three Goldstone modes only one survives 
reflecting the unchanged rotational symmetry around the 
chain axis. 

The spin- wave gaps at q = amount to i?o,i = 
198 GHz, E 0:2 = 295 GHz and E 0t3 = 410 GHz for 
CsMnBr 3 , which compares favorably with the experimen- 
tal values jl6],[l7]]. Note that this calculation has no free 
parameter to fit, since the dipolar energy is determined 
by the lattice constants. We also calculated the spin- 
wave gaps for RbMnBr3; neglecting crystal distortions 
we also obtain good agreement with the experiment pSj . 
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